Cosmology with a decaying vacuum 

K. Urbanowski^* and M. Szydlowski^* 

* University ofZielona Gora, Institute of Physics, ul. Prof. Z. Szafrana 4a, 65-516 Zielona Gora, Poland 
^ Jagiellonian University, Astronomical Observatory, ul. Orla 171, 30-244 Krakow, Poland 



cn 

o 

< 
O 



cr 



> 

O 



X 



Abstract. Properties of unstable false vacuum states are analyzed from the point of view of the quantum theory of unstable 
states. Some of false vacuum states survive up to times when their survival probability has a non-exponential form. At times 
much latter than the transition time, when contributions to the survival probability of its exponential and non-exponential 
parts are comparable, the survival probability as a function of time t has an inverse power-like form. We show that at this time 
region the instantaneous energy of the false vacuum states tends to the energy of the true vacuum state as 1/f^ for f ^ oo. 
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INTRODUCTION 

A discussion of the decay of the false vacuum was be- 
gan in two pioneer papers by Coleman and, Callan and 
Coleman [1, 2]. "The fate of the false vacuum" was dis- 
cussed there, namely the unstability of a physical system 
in a state which is not an absolute minimum of its en- 
ergy density, and which is separated from the minimum 
by an effective potential barrier It was shown, in those 
papers, that even if the state of the early Universe is too 
cold to activate a "thermal" transition (via thermal fluc- 
tuations) to the lowest energy (i.e. "true vacuum") state, 
a quantum decay from the false vacuum to the true vac- 
uum may still be possible through a barrier penetration 
via macroscopic quantum tunneling. Not long ago, the 
decay of the false vacuum state in a cosmological con- 
text has attracted renewed interest, especially in view of 
its possible relevance in the process of tunneling among 
the many vacuum states of the string landscape (a set of 
vacua in the low energy approximation of string theory). 
Since the work of Khalfin [3] it is known that for 
long times compared to the characteristic decay time of 
an unstable state (when the decay law has an exponen- 
tial form), the survival probability of such states is no 
longer described by an exponential function of time f 
but it decreases as f ^ oo more slowly than any expo- 
nential function of f . Krauss and Dent analyzing a false 
vacuum decay [4, 5] pointed out that in eternal inflation, 
even though regions of false vacua by assumption should 
decay exponentially, gravitational effects force space in 
a region that has not decayed yet to grow exponentially 
fast. This effect causes that many false vacuum regions 
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can survive up to the times much later than times when 
the exponential decay law holds. In the mentioned pa- 
per by Krauss and Dent the attention was focused on 
the possible behavior of the unstable false vacuum at 
very late times, where deviations from the exponential 
decay law become to be dominat. In our paper the atten- 
tion will be focussed on properties of decaying false vac- 
uum states from the point of view of the quantum theory 
of unstable states evolving in time and decaying. It will 
be shown that at late times the instantaneous energy of 
the false vacuum states tends to the energy of the true 
vacuum state as 1 /f ^ for f ^ oo. This means that in the 
case of such false vacuum states the cosmological con- 
stant becomes time depending. Next effective cosmolog- 
ical model with obtained decaying Lambda term will be 
studied in details. 



UNSTABLE STATES IN SHORT 

If \M) is an initial unstable state then the survival prob- 
ability, ^(f), equals ^(f) = |A(f)p, where A(f) is the 
survival amplitude, A (f) = (M|M;f), and fl(0) = 1, and 
[\M;t) = exp[-/f//] \M), H is the total Hamiltonian of 
the system under considerations. The spectrum, C7(//), 
of H is assumed to be bounded from below, (y{H) = 
[Emm,°°) and £■,„,„ > -oo. 

From basic principles of quantum theory it is known 
that the amplitude A(f), and thus the decay law ^M{t) 
of the unstable state \M) , are completely determined by 
the density of the energy distribution function (o{E) for 
the system in this state 



A(f) = 



Spec.(H) 



CO{E) e' 



-Et 



dE. 



(1) 



where (o{E) > and (o{E) —Q for E < Emm- From this 
last condition and from the Paley- Wiener Theorem it fol- 



lows that there must be (see [3]) |A(f) I > Ai exp[-A2f''], 
for |f| ^>oo. HereAi >0,A2 >0 and <q <1. 

This means that the decay law ^M(t) of unstable 
states decaying in the vacuum can not be described by an 
exponential function of time f if time f is suitably long, 
f — 7- oo, and that for these lengths of time ^M{t) tends to 
zero as f — >^ °o more slowly than any exponential function 
off. 

The analysis of the models of the decay processes 
shows that >^m(0 — exp[ — ^], (where Fm is the decay 
rate of the state \M)), to an very high accuracy at the 
canonical decay times t: From t suitably later than the 
initial instant to up to f 3> Tm == fi/^M and smaller than 
t — T, where T is the crossover time and denotes the time 
f for which the non-exponential deviations of A(f ) begin 
to dominate. 

In general, in the case of quasi-stationary (metastable) 
states it is convenient to express A(f) in the following 
form 

A(f)=A«p(f)+A„o„(f), (2) 

where A«p(f) is the exponential part of A(f), that is 
Aexp{t) = exp[-/f (^M - 5 Tm)], (Em is the energy of the 
system in the state \M) measured at the canonical de- 
cay times, A^ is the normalization constant), and A„o„(f) 
is the non-exponential part of A(f). For times t ^ Tm'- 

|A„p(f)l > l^«o«(f)l' 
The crossover time T can be found by solving the 

following equation. 



|A«p(f)|2 = |A„«„(f)|^ 



(3) 



The amplitude A„,,„{t) exhibits inverse power-law be- 
havior at the late time region: t ^T. Indeed, the integral 
representation (1) of A(f) means that A(f) is the Fourier 
transform of the energy distribution function co{E). Us- 
ing this fact we can find asymptotic form of A(f) for 
t —>■ °°. Results are rigorous. 

Let us consider co{E) having universal and general 
form. Namely let co{E)he of the form 



co{E) = iE-E,„i„)^ri{E) e Li(- 



(4) 



where < A < 1 and it is assumed that 77 {E,„i„) > and 
derivatives 77 '^^'(/i), (fc = 0, 1, . . . ,n), exist and they are 
continuous in [£■„,„, 00), and hmits lim£^£„_.„+ 77 (*)(£■) 



exist, lim^^oo {E — E„ 
tioned k, then 



■ r\W{E)=Q for all above men- 



A{t) 
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r(A + 1)770 
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(5) 



r(A+2)77^ 
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ENERGY AND DECAY RATE 

The amplitude A (f) contains information about the decay 
law ,'3^M{t) of the state \M), that is about the decay rate 
Ym of this state, as well as the energy Em of the system 
in this state. This information can be extracted from A(f ). 
Indeed if \M) is an unstable (a quasi-stationary) state 
then A (f) = exp[^(£'Af ^ ^Tm)?]. So, there is 



Em 



-M - 



ih 



dA{t) 1 



(6) 



in the case of quasi-stationary states. 

Taking the above into account one can define the "ef- 
fective Hamiltonian", /zm, for the one-dimensional sub- 
space of states Jif\\ spanned by the normalized vector \M) 
as follows 



def .^^ dA{t) 1 

dt A{t)' 



Hm — ih 



(7) 
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In general, Hm can depend on time f, Hm = hM{t)- One 
meets this effective Hamiltonian when one starts with 
the Schrodinger Equation for the total state space ,^ 
and looks for the rigorous evolution equation for the dis- 
tinguished subspace of states ,^u C ^. Details can be 
found in [6] and in [7]. Thus, one finds the following ex- 
pressions for the energy and the decay rate of the sys- 
tem in the state \M) under considerations, to be more 
precise for the instantaneous energy <§M{t) and the in- 
stantaneous decay rate, yuit), S'm = <^(f ) = 51 {hM{t)), 
jM = 7M(f) = -23 {hM{t)), where 5R(z) and 3 (z) de- 
note the real and imaginary parts of z respectively. 

Using (7) one can find that (^(f ^ Tm) — Em y^ <^(0), 
YM{t ^ Tm) — Fm, i. e., there is (^M{t) = Em at the 
canonical decay time. 

Starting from the asymptotic expression (5) for A(f) 
and using (7) after some algebra one finds for times 
t > T that 

/ZM(f )U=o ^ E^in + {-j)ci + {-jfc2 + ..., (8) 

where Ci ~ c* , / = 1,2,...; (coefficients c, depend on 
Co{E)). This last relation means that 

Mt) - £™„ + ^..., (forf»r). (9) 

(10) 

(Here 02 > 0). These properties take place for all unstable 
states which survived up to times t ^ T. From (9) it 
follows that lim(^oo S'M{t) = Emm- 



COSMOLOGICAL APPLICATIONS 

Krauss and Dent in their paper [4] mentioned earlier 
made a hypothesis that some false vacuum regions do 



survive well up to the time T or later. Let \M) — lO)-^"'*^, 
be a false,|0)" 



a true, vacuum states and E, 



false 




be 



the energy of a state corresponding to the false vacuum 
measured at the canonical decay time and E'™'' be the 
energy of true vacuum (i.e. the true ground state of 
the system). As it is seen from the results presented in 
previous Section, the problem is that the energy of those 
false vacuum regions which survived up to T and much 
later differs from E^ '*'', (see [8] and references one can 
find therein). 

Going from quantum mechanics to quantum field the- 
ory one should take into account among others a volume 
factors so that survival probabilities per unit volume per 
unit time should be considered. The standard false vac- 
uum decay calculations shows that the same volume fac- 
tors should appear in both early and late time decay rate 
estimates (see [4] ). This means that the calculations of 
cross-over time T can be applied to survival probabilities 
per unit volume. For the same reasons within the quan- 
tum field theory the quantity ^m [t ) can be replaced by the 
energy per unit volume Pm because these volume factors 
appear in the numerator and denominator of the formula 
(10) for hM{t)- Now, if one assumes that E'f^"" = E,„i„ 
then one has for the energy of the false vacuum state that 
at f > r. 



^false 



C2 



sr\t)c^E'r+'^... ^E, 



^ fahe 

'o 



(11) 



This property of the false vacuum states means that 

7true „„ . . „ ^2) 



Sf"''" {t) ^ E'q"" as f^oo. 



The basic physical factor forcing the wave function 
\M;t) and thus the amplitude A{t) to exhibit inverse 
power law behavior at f :^ T is a boundedness from 
below of o{H). This means that if this condition takes 
place and J^^ (o{E)dE < o°, then all properties of A(f ), 
including a form of the time-dependence at t ^ T, are 
the mathematical consequence of them both. The same 
applies by (7) to properties of hM{t) and concerns the 
asymptotic form of liMit) and thus of S'M{t) and 7M(f) at 

f >r. 

Note that properties of A(f ) and liMit) discussed above 
do not take place when g{H) = {—°°, +°°). The late time 
behavior of the energy of the system in the false vacuum 
state. 



jpfalse /■ \ 



Tifalse/^-^ ^ jj,true , ^2 
false 



-4..., for f>r, 



or, of the energy density Pq ' (f): 



false/^\ ^ ^true . "2 



pr'{t)c^pr+-^---, for f»r, 



(13) 



(14) 



(where di — d* , i —,l,2,.. .), is the pure quantum effect 
following from the basic principles of the quantum the- 
ory. The standard relation is p^"" = j^ = Pbare- 



FALSE VACUUM COSMOLOGY 

Let us take into account properties of the false vacuum 
described by the last relations and let us consider the 
cosmological model which is homogeneous and isotropic 
(possesses the Robertson - Walker metric) in which en- 
ergy density of dark sector is given by the relation which 
consists of the bare cosmological constant A and the term 
depending on the cosmological time t 



P =Pde=A- 



a^ 



f2 • 



(15) 



As it is well know the Friedmann-Robertson- Walker 
model has a first integral which does not depend on the 
pressure p and is given in the form called the Friedmann 
equation 



3//^ = Pde+Pn 



(16) 



dlna 



where H = ^^ is the Hubble parameter, t is the cosmo- 
logical time and a = fl(f ) is a scale factor. 

For simplicity we assume that space is flat and matter 
is in the form of dust (p ~ 0), so p = p,„ = Pm.o a^ 
describes a dependence of the energy density of the 
matter on a scale factor or a redshift z: 1 +z = a" 
Therefore (15) and (16) give 



-,-1 



PmflCl 
3//2 



-3 



A 
3H^ 



3//2 



= 1. 



(17) 



Next we introduce the state variables (x, y, z) defined as 
,/o — r = 



_ , / Pm.O" 



/Q.A == y = 



and 



z = y '2 ■ The quantities x, y, z has the cosmolog- 
ical sense of density parameters for a matter, the cos- 
mological constant and for the which mimics time de- 
pendence of the false decaying vacuum. The cosmolog- 
ical constant A is assumed to be positive. Above vari- 



ables, of course, satisfy the constrain x +>■ 



,2_ 



1. Us- 



ing this constrain after some algebra we obtain for state 
variables (x,y,z) dynamical system in the phase space 
(xjjjz) : x^ -\-y^ +zi^ ~ 1, (i.e. on the unit sphere), having 
the following form: 



dx 
dx 

dy 
dx 

dz 
dx 




(18) 
(19) 
(20) 



where z = ±-\/l— x^— y^ ^nd the cosmological time t is 
reparametrized to the new time parameter x which is a 
monotonic function off, i.e., x = Ina. The above system 



has three 2D invariant submanifolds: {x — 0}, {y — 0}, 
{z = 0} and the following critical points: 

1. x= 1,37 = 0, z = (Einstein-de Sitter universe), 

2. x = 0,y = I, z = (de Sitter stationary universe), 

3. z^ I, x = Q, y = (Milne universe), 
5V4-3a2, z = 5^302 (Einstein-de 



4. y = 0, x= iV4-3a2, z 

Sitter dominated by the time-dependent part of the 
cosmological constant). 



If we eliminate one state variables from the constraint 
condition then model dynamics constitute the 2D dynam- 
ical system and then right-hand sides of the system (18) 
— (20) are not given in the polynomial forms. In neigh- 
borhood of initial singularity the generic solution is rep- 
resenting ACDM universe dominating. There is a new 
critical point representing by a node and saddle point 
representing the Milne universe (without the horizon). In 
this case a — uqI and we have an unstable invariant sub- 
manifold, which is 2D if a ^ 1, and an unstable node if 
0: <C 1 . For the late time (f — > 00) the trajectories of the 
system are going to the Sitter universe like for the FRW 
model with the bare cosmological constant. 

From dynamical analysis of the system one can ob- 
serve how changes the critical points as well as its type. 
There are present bifurcation as the value of a parameter 
changes. For example the bifurcation value of ttbif for the 
critical point (4) is: ttbif — -??. 





Figure 2. A phase portrait of the system (18)-(20) for a = 
2. On the phase portrait there is one critical point at finite 
domain of saddle type and four critical points located at the 
boundary Sl„i + fi/v = 1 . 



CONCLUSIONS 

We have shown that at late times, t ^ T, the energy of 
the false vacuum behaves as follows: 

4""'{t)^E'o'-"^ + '^..., for f»r, (21) 
(c2 > 0), and similarly the energy density p^ '"'(f)- 

for f>r. (22) 



pf"'-it)^pr+^ 



Using the above property of the false vacuum we 
have studied the dynamics of the flat cosmological FRW 
model with decaying false vacuum parameterized by the 
cosmological time. 
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